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Parquet diagrams sum self-consistently Feynman graphs for the vertex function with 
all two-particle multiple scatterings. We show how the complete parquet equations for the 
Hubbard-like models can be integrated to a generating functional from which all thermo- 
dynamic quantities are derived via (functional) derivatives. An explicit Luttinger-Ward 
functional $[G; A, /C;i7] is constructed containing only the renormalized one-particle, G, 
irreducible, A, and reducible, /C, two-particle propagators as independent variational 
functions. The parquet approximation is proven to be a thermodynamically consistent, 
^-derivable theory obeying the necessary conservation laws. 



One of the major problems in the theory correlated electrons is to construct in a systematic and con- 
trolled way a consistent approximation interpolating reliably between weak- and strong-coupling regimes. 
The two extreme limits of weak and strong couplings in the archetypal Hubbard model can be described 
relatively well. The weak-coupling regime is governed by a Hartree-Fock mean field with dynamical fluc- 
tuations covered by Fermi-liquid theory. Extended systems at low temperatures are Pauli paramagnets 
with smeared out local magnetic moments. For bipartite lattices antiferromagnetic long-range order sets 
in at half filling and zero temperature at arbitrarily small interaction. In the strong-coupling regime the 
Hubbard model at half filling maps onto a Heisenberg antiferromagnet with pronounced local magnetic 
moments and the Curie- Weiss law for the staggered magnetic susceptibility, at least at the mean-field 
level. The spectral structure is dominated by separated lower and upper Hubbard bands and the strongly 
correlated system seems to be insulating even in the paramagnetic phase. 

However, it is the intermediate coupling, where the effective Coulomb repulsion is comparable with 
the kinetic energy and hence neither very weak nor very strong, that is of great interest for the the- 
orists as well as for the experimentalists. At intermediate coupling dynamical fluctuations control the 
low-temperature physics of interacting electrons and neither weak-coupling nor atomic-like perturbation 
theories are adequate. In this nonperturbative regime a singularity in a generic two-particle function is 
approached and we expect breakdown of the Fermi-liquid regime and a transition to an ordered state or 
eventually to a Mott insulator. 

Unfortunately there are only a few techniques for a quantitative analysis of the transition region 
between weak and strong coupling. Exact methods such as the Bethe ansatz or the numerical renormal- 
ization group can be applied only to one-dimensional or single- impurity models [^j^ . Numerical quantum 
Monte Carlo is good for thermodynamic properties at relatively high temperatures and is restricted to 
small samples ||^. Analytic methods are mostly of effective or mean-field nature Q. Systematic dia- 
grammatic expansions usually do not go beyond the fluctuation-exchange approximation (FLEX) that 
is known to fail at intermediate coupling Q. Although some improvements upon or alternatives to the 
FLEX approximation have been proposed, the vertex renormalization in these theories remain static 7|. 
The full dynamic and self-consistent vertex renormalization is contained first in the parquet diagrams [3|. 

Parquet diagrams were introduced to describe interaction of mesons more effectively Since then 
a number of attempts have been made to utilize the nontrivial renormalization scheme of the parquet 
algebra in condensed matter. Kondo effect ||T^, x-ray edge problem ||ll[], formation of the local magnetic 
moment fl^ are among the most well known applications. Inability to solve the parquet equations 
effectively has impeded broader application of the method. 

Parquet diagrams represent a systematic way of summation and renormalization of Feynman graphs. 
Instead of concentrating on the one-particle irreducible diagrams and the Dyson equation, the parquet 
approach sums two-particle diagrams contributing to vertex functions for which Bethe-Salpeter equations 
are constructed. The resulting algebra is much more complicated than that of the one-particle approx- 
imations. The two-particle Green functions obtained from the parquet approximation are suitable for 
spectral properties of the system, but it is cumbersome to gain thermodynamic properties out of them. It 
has been hitherto unclear whether the parquet approximation forms a thermodynamically consistent 
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derivable approximation fulfilling the necessary conservation laws . Only if we construct a generating 
Luttinger-Ward functional in closed form from which all the thermodynamic quantities can be derived 
via derivatives with respect to auxiliary sources, we can be sure thermodynamic relations and sum rules 
are fulfilled. 

The author proposed recently a simplification of the full parquet approximation by summing only 
two singlet two-particle channels yielding most divergent diagrams in the critical region of the metal- 
insulator transition p4| . A generating functional for such an approximation was derived via the linked- 
cluster theorem. The question is whether also the unabridged parquet equations can be integrated to 
a generating thermodynamic potential. It is the aim of this paper to show how the Luttinger-Ward 
functional and the grand potential can be constructed for the complete parquet equations, i. e. with 
all two-particle irreducible diagrams, in lattice models of interacting electrons with a local (Hubbard) 
interaction. 

There are three topologically inequivalent definitions of a two-particle irreducibility. It may be defined 
according to cutting electron-hole or electron-electron pair propagation, or according to cutting polariza- 
tion bubbles shielding the electron-electron interaction. Each possibility defines a two-particle channel 
of multiple pair scatterings characterized by a different binding of independent variables in the vertex 
functions. A general two-particle quantity will be denoted in this paper as in Fig. 1. Each two-particle 
function carries three independent four-momenta and two spin indices. We generally denote the fcrmionic 
four-momenta k = (k, iujn) and the transferred bosonic ones as q = (q, iiym), where cun — (2f^ + 1)ttT and 
i^m = 2'rmrT are the respective Matsubara frequencies at temperature T = l3~^. 

It is convenient to introduce a matrix notation in the spin indices to distinguish different two-particle 
channels in the parquet diagrams. We define a 2 x 2 matrix for the generic two-particle function Xa-a'- 



X = ( ^TT 



(1) 



We speak about singlet and triplet contributions to a two-particle function if the spins of the involved 
fermions are antiparallel or parallel, respectively. 

We define three matrix multiplication schemes for two-particle quantities 



X •¥ 



X oY 



X*Y 



, (fc, k'- q)^^J2 ^-'(fc, k'; q")Y,,, (fc + q" , k' q" ■ q - q") , 

9" 

, ik,k';q) = -^^X„„'{k,k' + q":q - q")Y„„'{k + q - q" ,k';q") , 

q" 

{k,k'-q) = ^ J2 X,,„{k,k";q)Ya"a'{k\k'-q) 

a"k" 



(2a) 
(2b) 
(2c) 



representing summations in the three inequivalent two-particle channels, electron-hole, (e/i), electron- 
electron, (ee), and interaction, ([/), respectively. We see that the variables of the two-particle functions 
are convoluted differently in inequivalent channels. Note that only the interaction channel mixes the 
singlet and triplet contributions. 

We decompose the full two-particle Green function into a sum of always reducible and irreducible 
projections onto each inequivalent channel 



tCaa' {k, k';q) ^ IC"^, (fc, k';q)+ I^^, (k, k'; q) 



(3) 



where a = eh, ee, U refers to a two-particle channel. 

The parquet diagrams can at best be represented graphically. Having in mind the above introduced 
notation and a general convention that double primed variables are summed over, we can write 
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ak 



ak + q 



k + q" 





a'k' + q 



ak + q 



k' + q" 



k + q- q" 



(4a) 





cr'k' + q 



k' + q" 



(4b) 



ak 



ak + q 



KF , 



k" 



( 



a'k' 



a'k' + q 



) 



k" + q 



(4c) 



We labeled only the intermediate four-momenta on the right-hand side of the parquet equations since the 
endpoints are the same as those on the left-hand side. 

Equations (|^) are the Bethe-Salpeter equations and constitute one set of the full parquet algebra. 
To complete it we must add relations connecting the two-particle reducible, /C", and irreducible with 
two-particle bubbles, A", functions. To this purpose we introduced horizontal {eh and ee) and vertical 
([/) two-particle bubbles 



G^'j^fc, k'; q) = G,{k + q)G,. {k' + q) , G^^>r{k, k'; q) - G,, (fc')G,, (fc' + q) 



{2)v, 



(5) 



The functions A" can then be defined as [11 5| 



K.'ik,k';q) 



C/(5,,,_,-|-AW(fc,fc';g)+ ^ lcCik,k';q) 



Gi'j"(fc,fc';g) 



(6) 



where is a sum of all Feynman diagrams simultaneously irreducible in each two-particle channel. 

The sum contains only higher-order diagrams where three and more particles are multiply interconnected. 
In the usual treatment with only two-particle multiple scatterings this irreducible part is neglected. We 
hence put A/^^' — 0. 
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The parquet algebra for the two-particle function is now complete. Equations and are to be 
integrated to a generating functional. We first construct two functionals being integrals of equations (Q) 
and (|^). We can exclude the reducible parts /C" from (Q) and the irreducible ones A" from (|). If we 
denote the generating functionals for (H) and (0) <f>A and respectively, we demand 



A" 



K 



(7) 



where A'J^, = X^,^ and X'^^, = X'^^, is a matrix transposition in the vertical and horizontal channels, 
respectively. The functional <I>a will be a sum of functionals generating each of the parquet equations 
in (^), i. e. each channel will have its own generating functional. Recalling the matrix notation and the 
multiplication rules (H) we can integrate each parquet equation to the respective generating functional 



ak a' k' 



2/32A/-2 



(fc,fc';0) 



(8a) 



2/3W' 



EE 

ak cr'k' 



2 



In 1 + A"^ 



(8b) 



2/32 A/-^ 



■EE 

<T kq 



A^-^A'^^A^ _ln(i + A^* 



1- 



'U 



(l + A^*) 



{k,k]q) 



(8c) 



The symbols A», Ao, A* in the logarithms set the appropriate multiplication for the power-series definition 
of the matrix function. 

It is much easier to integrate equation dfl) . It is straightforward to verify that 



2/32 A/-: 



ak a' k' 



2/32 AA^ 



■EE [{yS + K."''' + K.''^'^ G'^'^^^^^G^^)"] {k,k,q), 

a kq 



(9) 



where Jo-o-' = <^(t',-o-- Note that when we multiply functions from the electron-hole and electron-electron 
channels the multiplication rule can be taken from either channel without influencing the result. Con- 
necting the incoming and outgoing fermion lines in the product is uniquely determined by the charge 
conservation. 

The Luttinger-Ward functional must finally be stationary with respect to first variations of the two- 
particle functions. This can be achieved by subtracting the mixed trA/C term. Moreover, the functional 
must generate the self-energy via a variation with respect to the full one-electron propagator. We hence 
must add a purely one-electron part to the generating functional. The one-electron contribution is just a 
second-order term in the free-energy expansion. The parquet diagrams sum multiple scatterings of pairs 
of quasiparticles starting with the next term beyond second order. We choose the correct overall sign of 
the generating functional from the self-energy, e. g. in second order. The Luttinger-Ward functional for 
the parquet diagrams finally reads 



<i>[G;A,/C;C/] = -1^5:5:[A-.£-G(^)^; 

a kq 



{k,k,q) 



+ 



2/32 AT' 



■EE 



-luSG^'^^'' • U6 + 1^^ • K"'' -f A"" o ^-^^ I g(2)'' 



(fc,fc';0) 



-<^K [G; /C; U] - $f [A-^''] - [A'^'^] $^ [A^] 



(10) 



where we introduced in each functional its explicit dependence on the variational one- and two-particle 
functions and the bare interaction U . 
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The grand potential generating the complete thermodynamics of the parquet approximation is con- 
structed from the Luttinger-Ward functional and a free-electron term with the Hartrec contribution . 
We hence add also the particle densities as additional variational parameters and obtain for a fixed 
^„ = -\- ah where is a chemical potential and h an external magnetic field 

^Vl[n^,ni] S, G; A,/C] = X! hi [icj„ + ^i^ - e(k) - Un^„ - Scr(k, iLj„)] 

+Ga{KiuJn)^a{\^,iuJn)]-Un^ni + '^[G-A,lC-U] . (11) 
Here n|, n^; E, G; A, K, are independent variables and complex functions. Their physical values are cho- 



sen from stationarity of the grand potential (11) with respect to variations of its independent vari- 
ables/functions. There are three pairs of Legendre conjugate variational variables in (jl^). The Hartree 
parameters n| and n|, the one-electron functions Ticik) and G„{k), and finally the two-particle vertex 
functions A"^, (A:, k'\ q) and IC"^,{k, k';q). 

The stationarity equations for the Hartree parameters lead to a one-electron sum rule 

1 v-^ gii^TiO+ 

~ AA^ ^ iojn + fia- e(k) - Un-a - Sa(k, zw„) ' ^^^^ 

The stationarity equations for the one-particle functions, i. e. the one-electron propagator and the self- 
energy, measuring in this formulation corrections to the Hartree approximation, read 

Gcr(k, iw„) = [iw„ + - e(k) - C/n_cr - S^(k, iwn)]^"^ , (13a) 
^-(^) = E E ^-(^ + «)'^--' + 1^ ~Q)Ga'{k')Ga'{k' + q) . (13b) 

ct' k' q 

Finally, the stationarity equations for the two-particle functions are just the parquet equations (jj) and 
(|^) with A/o-cr' = 0. The full two particle Green function is a sum of a second-order contribution and the 
reducible parts from each two-particle inequivalent channel 

/C,,. {k, k'-q)^ USa,,^a + K.t' {k, k'- q) + (fc, k'- q) + IC^,, {k, k' ; q) . (14) 

Thermodynamics of the parquet approximation is completely determined by (|^)-(|l4|) into which ap- 
propriate auxiliary external sources are introduced. All thermodynamic relations are hence fulfilled in 
the parquet summation scheme. Moreover due to the full self-consistency in the one- and two-particle 
propagators, the parquet approximation obeys momentum and energy conservation laws p^ ]. However, it 
is unclear whether the parquet approximation conserves also charge as the only generator of the electro- 
static energy. The electrostatic energy (due to the Coulomb interaction U) of the solutions must entirely 
be generated by the actual distribution of charge carriers (electron density n). This conservation law 
induces a Ward-like identity connecting a variation in the electrostatic energy with a variation in the 
mass distribution . The identity relates a thermodynamic correlation function defined as a derivative 
of the one-electron density with respect to an auxiliary field (electro-chemical potential) and the full 
two-particle Green function JC-^i from (p^. 

The correlation function describing the variation of the electrostatic energy can be derived from 
the generating functional if we introduce a small perturbation of the local Coulomb interaction U 
U -(- SU{q). We define the correlation function as 



S^[G;A,IC;U 
C{q) = ■ 



5U=0 



5U{q) 

= E fc'; 9) {^Ti {k, k'; -q) + K.U {K k'; q)} G^f {k, k'; -q) . (15) 

The static part of this dynamic correlation function is important for the thermodynamic stability of a 
particular solution (phase). This function must be negative in a stable phase 
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C(q,0)<0 (16) 

as a consequence of a repulsive character of the Coulomb interaction The correlation function C(q, 0) 
is from definition connected with the magnetic susceptibility x(q) and the compressibility ^(q) 

^(^'0) = ;;^E^^"*^'''"''^^''K"iT^ji) - («iT)(^ji)] ^ J [«(q) -x(q)] , (17) 



which together with (|16|) indicates that in tight-binding models with a local electrostatic interaction the 
spin fluctuations are stronger than the charge ones. 

Equality is rigorously fulfilled in an exact solution. The Ward identity due to conservation 
of sources of the electrostatic energy demands that sums of both sides of ( p^ over momenta must 
equal, otherwise there are spurious sources of the electrostatic potential in the system. Approximate 
theories with only one-electron renormalizations break equality ( p7| ) qualitatively, i. e. the left- and right- 
hand sides lead to incompatible phase diagrams and it is unclear how to define two-particle functions 
consistently. The only way to improve upon this discrepancy is to introduce dynamical vertex (two- 
particle) renormalizations. It is clear that the exact equality in ( p^ can be achieved in self-consistent 
theories only in a complete solution where arbitrarily large clusters of particles are involved. In order to 
produce unambiguous results for two-particle correlation functions it must be required from approximate 
theories that the left- and right-hand sides of (|l^) generate the same divergences and lead to equivalent 
phase diagrams. 

Although the parquet diagrams contain an advanced dynamical renormalization of the vertex func- 
tions, skeleton three-particle d iagr ams are not involved. Higher-order skeleton diagrams are generated by 



derivatives of the self-energy (13b) with respect to the electro-chemical potential when constructing the 
susceptibility or the compressibility. However, skeleton diagrams beyond the parquet approximation are 
fully renormalized in the two-particle sector, i. e. they use appropriate projections of the full two-particle 
Green function JC^a' instead of the bare interaction U. The lowest-order contribution to A/g-o-' is shown in 
Fig. 2. Unlike other simpler approximations the two-particle functions in the parquet approach can show 
only integrahle singularities. A singularity can arise in one or more two-particle channels simultaneously. 
Only one bosonic variable, the conserved transfer four-momentum in the multiple scattering events, is 
relevant in each channel. It is for the electron-hole, electron-electron, and interaction channel k — k' , 
k + k' + q, and q, respectively. The divergences in two-particle functions are smeared out (regularized) 
whenever one integrates over the variables in which the singularity arises. We realize that three-particle 
skeleton diagrams as in Fig. 2 contain integrals over all relevant two-particle four-momenta simultane- 
ously and are hence finite. Unless qualitatively new singularities arise in three-particle (or higher-order) 
skeleton diagrams, the parquet approximation is stable with respect to fluctuations and the left- and 
right-hand sides of (^^ lead to the same instabilities and hence to qualitatively the same phase dia- 
gram. Appearance of new singularities in three-particle and higher-order skeleton diagrams, not seen at 
the two-particle level, would indicate in the renormalization-group language that higher-order effective 
interactions are relevant. Neither the Bethe ansatz nor the renormalization group suggest relevance of 
interactions beyond the two-particle ones. 

Parquet diagrams are a systematic approximation summing all two-particle irreducible Feynman 
graphs with fully renormalized one-electron propagators. The idea behind this summation scheme is 
to replace an expansion in the bare interaction strength U by skeleton diagrams with renormalized two- 
particle functions. This is a realization of the field-theoretic renormalization of perturbation theory 
close to critical points. Parquet approximation is the simplest theory with two explicit vertex functions 
containing minimally one-loop renormalization of the critical behavior. However, one cannot reduce the 
singular two-particle Green function to a static running coupling constant as in classical phase transitions. 
As well as we cannot reduce renormalized perturbation theory to a weak-coupling loop expansion. Due 
to a complicated structure of dynamical fluctuations at quantum phase transitions in itinerant systems, 
one has to keep dynamical two-particle Green functions in the parquet equations and their extensions. 
Parquet diagrams offer a comprehensive basis for a reliable assessment of the critical behavior at quantum 
phase transitions with singularities in two-particle Green functions. 

Extensions of the parquet approximation contain only skeleton diagrams contributing to Al^a-' free 
of one- and two-particle insertions. The interaction strength is everywhere replaced by the two-particle 
Green function. The parquet approximation offers an overall reliable solution whenever the norm of 
the full two-particle function || /CG^^^ \\<^ 1. Since the parquet equations allow only for integrable 
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singularities in the two-particle Green functions, their solution may obey the reliability criterion even at 
intermediate and strong coupling. At least || ICG^'^^ \\< 1 is guaranteed in the parquet approximation. 
This feature makes it an appropriate tool for studying the transition region between the weak-coupling, 
Fermi-liquid solution and the strong-coupling, local-moment phase in the lattice models with an effective 
local (Hubbard) interaction. 

To conclude, we integrated the parquet diagrams for two-particle Green functions to a generating 
Luttinger-Ward functional. We demonstrated that the parquet summation scheme belongs to thermo- 
dynamically consistent and conserving ^-derivable approximations in Baym's sense if the two particle 
functions are generated from appropriate functional derivatives of the self-energy. The difference be- 
tween the thermodynamic (variational derivatives) and diagrammatic definitions of two-particle functions 
(Aslamazov-Larkin diagrams |l^), reflecting deviations from exact identity (|l^), are expected to vanish 
in leading asymptotic order at two-particle critical points. 

Although the explicit grand potential does not help to resolve the complicated algebra of dynamical 
variables coupled in the parquet equations, it enables one to control systematic expansions beyond and 
eventual simplifications within the parquet approximation. The parquet diagrams and their derivatives 
have been presented in a new light from which the role, relevance, and systematics of the parquet-type 
graphs have become more transparent. 
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ak + q 



a'k' 



a'k' + q 



FIG. 1. Generic two-particle function with three independent four-momenta and a defined order of incoming 
and outgoing fermions. 



k + q + q" 



ak + q 




oc k" 



a'k' + q 




k' + q" 



FIG. 2. Lowest-order contribution not included in the parquet approximation. The bare interaction is replaced 
by the full two-paxticle function from the parquet approximation. Double primed variables are summation indices. 
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